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ABSTRACT 



Precise measurements were conducted in continuous flow seawater mesocosms 
located in full sunlight that compared metabolic response of coral, coral-macroalgae 
and macroalgae systems over a diurnal cycle. Irradiance controlled net photosynthe- 
sis (P n et)> which in turn drove net calcification (G ne t), and altered pH. P net exerted 
the dominant control on [COj - ] and aragonite saturation state (f2 arag ) over the diel 
cycle. Dark calcification rate decreased after sunset, reaching zero near midnight 
followed by an increasing rate that peaked at 03:00 h. Changes in f2 arag and pH lagged 
behind G ne t throughout the daily cycle by two or more hours. The flux rate Pnet was 
the primary driver of calcification. Daytime coral metabolism rapidly removes dis- 
solved inorganic carbon (DIC) from the bulk seawater and photosynthesis provides 
the energy that drives G ne t while increasing the bulk water pH. These relationships 
result in a correlation between G net and f2 arag , with £"2 arag as the dependent variable. 
High rates of H + efflux continued for several hours following mid-day peak G ne t 
suggesting that corals have difficulty in shedding waste protons as described by the 
Proton Flux Hypothesis. DIC flux (uptake) followed P ne t and G ne t and dropped 
off rapidly following peak P ne t and peak G n et indicating that corals can cope more 
effectively with the problem of limited DIC supply compared to the problem of 
eliminating H + . Over a 24 h period the plot of total alkalinity (Aj) versus DIC as 
well as the plot of G net versus £2 arag revealed a circular hysteresis pattern over the diel 
cycle in the coral and coral-algae mesocosms, but not the macroalgae mesocosm. 
Presence of macroalgae did not change G ne t of the corals, but altered the relation- 
ship between £2 arag and G ne t- Predictive models of how future global changes will 
effect coral growth that are based on oceanic £2 arag must include the influence of 
future localized P net on G ne t and changes in rate of reef carbonate dissolution. The 
correlation between f2 arag and G ne t over the diel cycle is simply the response of the 
C02-carbonate system to increased pH as photosynthesis shifts the equilibria and 
increases the [CO3 - ] relative to the other DIC components of [HCO^~] and [CO2L 
Therefore f2 arag closely tracked pH as an effect of changes in P ne t, which also drove 
changes in G ne t- Measurements of DIC flux and H + flux are far more useful than 
concentrations in describing coral metabolism dynamics. Coral reefs are systems that 
exist in constant disequilibrium with the water column. 
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INTRODUCTION 

Recent field experiments have identified the need for accurate metabolic measurements 
on coral reefs at short time intervals in order to detect subtle aspects such as phase lags 
between the concentrations and flux rates of major metabolic parameters over the diurnal 
cycle (e.g., Shamberger et ah, 2011; McMahon et ah, 2013). There is a need to test the 
assumption that night calcification is very low and constant and that calcification is 
limited by inorganic carbon concentration in the form of CO 2- . Finely detailed diurnal 
metabolic measurements are difficult to obtain in field investigations. For example, Falter 
et ah (2012) based their extensive study at Nigaloo Reef, NW Australia on a total of 13 
summer measurements of calcification taken at various times over 12 days between 
08:00 and 18:00 with one night measurement at 21:00. The following winter a total of 
1 1 more measurements were made over 4 days with no night measurements. These data 
are valuable, but do not provide the diurnal resolution that is needed to fully describe 
coral and coral reef metabolism. Thus we undertook a laboratory investigation that would 
provide such data. 

BACKGROUND INFORMATION 

The term £2 arag is defined as: 
[Ca 2+ ][C0 2 ~] 

q — - _ 3— m 

11 'arag — „ V 1 / 

where K sp is the solubility constant of aragonite. The [Ca 2+ ] in normal present-day oceanic 
seawater is essentially constant at 10.3 mmol kg -1 SW, normalized to salinity. Likewise, 
K s p is a constant (at a given temperature, pressure, and salinity), so in shallow oceanic 
waters f2 arag is directly proportional to [CO 2- ]. Changes in seawater pH shift the equilibria 
between the various forms of dissolved inorganic carbon (DIC) as follows: 

C0 2 + H 2 0<^H + + HCO~ <^2H + + C0 2 ~. (2) 

Calcification inevitably produces an excess of H + and thus reduces total alkalinity (At) by 
two moles for every mole of CaC03 precipitated (Kinsey, 1978; Smith & Kinsey, 1978). The 
correct equations for calcification are as follows: 

Ca 2+ + (C0 2 + H 2 0) & CaC0 3 + 2H+ (3) 

t 

Ca 2+ + (H + + HCC>3") <^CaC0 3 + 2H+ (4) 

$ 

Ca 2+ + (2H + + C0 2 ") <^CaC0 3 + 2H+ (5) 

Equations 3-5 are written in two dimensions with a red arrow showing the relationship 
between the carbonate species (in parentheses) that shift with the changes in [H + ] 
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described as Eq. (2 ). Dissolution is the reverse of the calcification reaction. Net calcification 
(G ne t) is the sum of calcification (positive flux) and dissolution (negative flux). When 
the equations are written correctly in this manner the importance of protons becomes 
apparent with two moles of H + produced for every mole of CaCC>3 precipitated regardless 
of which form of dissolved inorganic carbon (DIC) is involved. 

The following equations describe photosynthetic carbohydrate formation from the 
various available CO2 species: 

(C0 2 + H 2 0) & CH 2 0 + 0 2 (6) 

$ 

(H + + HCO-)^CH 2 0 + 0 2 (7) 

$ 

(2H + + C02-)^CH 2 0 + 0 2 . (8) 

The photosynthesis equations are also written in two dimensions with the red arrows 
showing changes in distribution of species that occurs (Eq. (2)) with shifts in pH. Note 
that photosynthesis increases pH (lowers [H + ]) while the reverse reaction of respiration 
decreases pH (increases [H + ]). Net photosynthesis (P ne t) is the sum of photosyn- 
thesis (positive flux) and respiration (negative flux). Unlike calcification-dissolution, 
photosynthesis-respiration does not alter total alkalinity (At). 

In sum, photosynthesis and calcification both lower the seawater DIC, while respiration 
and CaCC>3 dissolution raise DIC. Only the precipitation or dissolution of CaC03 
significantly alters Aj. Consequently, changes in [At] can be used to calculate calcification 
and dissolution rates (G ne t), and is widely used in this regard. Photosynthesis and respira- 
tion can radically alter [H + ] and thus can alter relative concentration of CO^.HCOJ" 
and CO2. Coral calcification is a biological process that is heavily influenced by the 
associated processes of photosynthesis and respiration (P ne t) that modify pH. Protons 
can be considered a waste product of calcification Eqs. (3)— (5) and O2 a waste product of 
photosynthesis Eqs. (6)-(8). 

METHODS AND MATERIALS 

A mesocosm experiment was undertaken in order to precisely measure the changes in 
bulk sea water chemistry and material flux caused by coral and algae metabolism over a 
diurnal cycle. The experiment was conducted in the flow-through mesocosm system at the 
Hawaii Institute of Marine Biology, Kaneohe Bay, Oahu, Hawaii. The mesocosm system 
has been described previously in detail (Jokiel et ah, 2008; Andersson et ah, 2009; Jokiel, 
Bohr & Rodger s, in press). The fiberglass mesocosm tanks were located in full sunlight 
and supplied with flowing seawater pumped from approximately 2 m depth at the edge 
of the Coconut Island coral reef. Each mesocosm received a flow of approximately 7.5 
to 8.5 1 min -1 resulting in a turnover rate of approximately 1 h. Solar input at the site 
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Table 1 Mesocosm volume and biomass of the coral Montipora capitata and the macroalga 
Gracillaria salicornia. 







Tank dimensions 


in cm 




Mesocosm 


Biomass 


Length 


Width 


Water depth 


Volume (1) 


Coral only 


7,555 g buoyant weight coral 


117 


117 


38 


520 


Coral plus Algae 


7,555 g buoyant weight coral, and 
3,151 g of Gracillaria salicornia 


117 


117 


38 


520 


Algae 


3,151 g of Gracillaria salicornia 


117 


117 


35 


479 



was monitored with a LiCor Brand Quantameter (Li-Cor Inc., Lincoln, NE, USA), which 
measured photosynthetically active radiance (PAR) between 400 nm and 700 nm. 

Three continuous flow mesocosms were used for this experiment. The first mesocosm 
("Coral only") was loaded with 7.6 kg buoyant weight of the reef coral Montipora capitata 
for close to 100% coverage of the bottom (Table 1). This buoyant weight translates into 
11.5 kg dry skeletal weight {Jokiel, Maragos & Franzisket, 1978). The second mesocosm 
("Coral plus Algae") contained the same weight of live coral plus 3.1 kg wet weight of the 
macroalgae Gracillaria salicornia. The third mesocosm ("Algae only") was loaded with the 
same weight of the macroalgae. A small biomass of calcifying organisms was present on 
the macroalgae in the form of epiphytes that were not removed. Dead skeletal material 
and sediment were excluded from all three mesocosms to reduce the complicating effect 
of decalcification of carbonates and related processes on G ne t (Murillo, Jokiel & Atkinson, 
2014). The organisms in the mesocosms were allowed to acclimate to mesocosm conditions 
for one week prior to the experiment. 

This experimental design allows the three treatments to be run simultaneously which 
eliminates the between-treatment variance due to solar irradiance, temperature and other 
factors that would occur if they were run at different times. The functioning of each 
community under identical conditions serves as a control and as a contrast to the other 
two communities. Response of these communities to the diurnal irradiance cycle and 
variation in other factors is non-linear, so the data were analyzed by graphical comparison 
and integrated response of 24 h P net and G ne t- 

The seawater inflow enters at the bottom in the center of each mesocosm (see Smith 
et ah, 1977 for details on the mesocosm system) at a vertical angle, which ensures a 
uniform and well-mixed system. Maintaining seawater inflow of reef water at constant 
rates insured that natural fluctuations in seawater chemistry observed on the adjacent reef 
during the diel cycle was preserved during the experiment. Water chemistry was sampled 
hourly from 06:00 on 24 April to 07:00 on 25 April 2012 at the inlet and outlet and the 
flow rate recorded. During the experiment the flow rate was precisely measured every 
hour for each mesocosm and inlet-outlet chemistry was determined. Temperature and 
salinity were measured with a YSI Brand Model 30 salinity-conductivity-temperature 
meter (±0.1 °C; ±0.1 ppt). Dissolved oxygen (DO) was measured with a YSI Brand 
Model 57 Dissolved Oxygen Meter (±0.2 mg 1 _1 ), and pHnbs with an Accumet AP72 
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pH/mV/temperature meter verified spectrophotometrically using m-cresol purple dye 
according to SOP 7 (Dickson, Sabine & Christian, 2007). Ax was measured using a Titrino 
Model 877 titrator system. Alkalinity samples were equilibrated to 25 °C and run within an 
hour of being taken. Accuracy and precision of the titrations was confirmed with certified 
reference materials (CRM Batch 129) from the Dickson Laboratory, Scripps Institution of 
Oceanography which verified that our measurements of Aj were accurate to within 0.18% 
of the Dickson CRM value. All carbonate parameters were calculated using the program 
C02SYS (Pierrot, Lewis & Wallace, 2006) and stoichiometric dissociation constants defined 
by Mehrbach et al. (1973) and refit by Dickson & Millero (1987) . 

The mesocosms are well mixed systems and can be represented by a simple box model 
Andersson et al. (2009). Changes in total alkalinity (A j) are attributed to calcification or 
carbonate dissolution Murillo, Jokiel & Atkinson (2014). The net calcification (G ne t) can be 
easily calculated from Eq. (9) by measuring the amount of material in the inflowing (Pin) 
and outflowing (P ou t) seawater, and the change per unit time in the mesocosm seawater 
between consecutive sampling times. For G ne t the resulting value is divided by 2 because 
two moles of Aj are produced for every mole of CaC03 precipitated (calcification) or 
removed (dissolution) as shown by Eqs. (3)-(5). 

F m A T -F out A T - d -£ 

^net — ■ (?) 

P net was measured in a similar manner using O2 concentration, with one mole of carbon 
being produced for every mole of O2 produced as shown by Eqs. (6)-(8). Likewise DIC and 
H + were calculated using Eqs. (11) and (12). 

Pnet - Pin0 2 - P ou t0 2 - — - (10) 

dt 

dDIC , 

DIC net = P in DIC - P out DIC - — — (11) 

dt 



dH+ 
dt 



= Pi n H + — P ou tH + — . (12) 



RESULTS AND DISCUSSION 

All of the measurements taken during this experiment are shown in Table 2. 

Calcification over the 24 h period (Fig. 1) shows the diurnal pattern related to 
irradiance, light enhanced calcification and dark calcification. Values for G net are high 
due to the large biomass of live coral, high solar irradiance in the shallow mesocosms and 
absence of sediment or dead carbonate skeleton which are subject to dissolution. G ne t 
in the "Corals only" and "Corals plus Algae" treatments track each other closely. Light 
saturation of calcification did not occur up to the maximum irradiance which exceeded 
1,500 (xmole photons m -2 s . This value is many times higher than that supplied by the 
artificial light typically used is most laboratory studies of coral calcification. Calcification 
rate is very low in the "Algae only" treatment due to low biomass of calcifying organisms, 
which are made up of various calcifying epiphytes. Low-levels of dark calcification occur at 
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Figure 1 Diurnal net calcification rate (G n et) and irradiance for the three mesocosms. 

night. There is a drop in calcification to zero around midnight with a dark calcification rate 
peak at approximately 03:00 h. 

The relationship between G ne t and 9. arag 

The linear regression of G net plotted as a function of f2 arag (Fig. 2) has become a 
widespread method of describing coral and coral reef calcification. A significant statistical 
relationship is obtained, with substantial variance that is generally assumed to be 
largely sampling error, or the result of other factors which influence G ne t- However, the 
explanation appears be more complex, as will be discussed in the following sections. 

Diel hysteresis, phase lags and night calcification patterns 

McMahon etal. (2013) quantified G ne t in a healthy coral reef lagoon in the Great Barrier 
Reef during different times of day. Their observations revealed a diel hysteresis pattern 
in the G ne t versus £2 arag relationship. This phenomenon can be demonstrated by labeling 
the points in Fig. 2 with the hour of day as shown in Fig. 3. The diel pattern moves from 
the lower left quadrant early in the day toward the upper right through mid-day and then 
back to the lower center during the night before returning to the lower left quadrant at 
first light. The pattern is nearly identical for the "Corals only" mesocosm (Fig. 3A) and the 
"Corals plus Algae" (Fig. 3B), which tracked each other closely (Fig. 1). The "Algae only" 
mesocosm did not show this pattern. The linear regression for the G ne t vs. f2 arag data for 
the mesocosm with coral (Fig. 2) accounted for part of the variance (R 2 = 0.40). A linear 
regression does not adequately describe the variance resulting from the diel pattern. 
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Figure 2 Net calcification rate (G ne t) plotted as a function of ft a rag f° r (A) "Coral only" mesocosm, 
(B) "Coral and Algae" mesocosm and (C) "Algae only" mesocosm. 



Cyronak et al. (2013) used chambers to measure in situ benthic solute fluxes at three 
different advection rates at Heron Island lagoon, Australia and observed a strong diurnal 
hysteresis pattern similar to that in Fig. 3. They suggested that diel hysteresis is caused by 
the diurnal interaction between photosynthesis and respiration. The data did not follow 
a trend consistent with the £2 arag of the water column being the main driver of carbonate 
precipitation and dissolution. Instead, carbonate precipitation and dissolution in these 
sediment communities is linearly correlated to the rates of photosynthesis and respiration 
CPnet) occurring over the same time period. 

Phase shifts 

Evaluation of phase relations for the various parameters listed in Table 2 can be facilitated 
by scaling each variable on a 0 to 1 scale. The normalized value of a; for variable A in the ith 
row was calculated using the equation: 

Normalized value (a;) — — mm (13) 

Amax A m ; n 

where A m ^ is the minimum value for variable A and A max is the maximum value for 
variable. Figure 4 summarizes the results for the variables most often considered in the 
literature (pH, fi arag , P n et and G net ). 
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Figure 3 Net calcification rate versus $2 a rag- Net calcification rate versus £2 ara g with each point labeled 
with hour of the day revealing the clockwise coral reef diel hysteresis pattern for: (A) "Coral only" 
mesocosm, (B) "Coral plus Algae" mesocosm and the (C) "Algae only" mesocosm, which showed no 
pattern. 



Figure 4 shows that peak pH and £"2 arag lag behind G net throughout the daily cycle by 
two or more hours. The figure also shows that peak G ne t follows P ne t during daylight 
photosynthetic hours with a reverse during the nighttime hours. Shamberger et al. (201 1 ) 
reported that £2 arag lags behind G ne t on the reefs of Kaneohe Bay, Hawaii. McMahon et 
al. (2013) reported that peak G ne t rates occurred 2-3 h before the £2 arag maximum on 
a healthy coral reef on the Great Barrier Reef. Thus £2 arag (along with closely correlated 
[CO3 - ], pH and [DIC]:[H + ] ratio) is not the primary driver of coral calcification over 
a diurnal cycle. The paradigm that f2 arag correlates with G ne t on a global scale must 
be tempered with the caveat that other processes have a much greater influence on 
calcification on smaller spatial and temporal scales. The data presented above show that 
diurnal irradiance drives P ne t> which in turn drives G ne t, and which alters pH, which 
controls [CO3 - ] and £2 arag as well as other variables based on concentration such as the 
ratio of [DIC] to [H + ]. A better understanding of this hierarchy will resolve many of the 
contradictions in the literature on coral reef calcification. 

Night calcification 

Laboratory studies show that coral calcification continues in darkness, but at a lower rate 
than observed in light enhanced calcification (Schneider & Erez, 2006). Night calcification 
rates have generally been assumed to be low and constant at night, although this assump- 
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06:00 12:00 18:00 00:00 06:00 

Time of Day 

Figure 4 Normalized pH, S2 a rag> Pnet and G n et values for the three mesocosms versus time of day using 
data from Table 2 and Eq. (9). Arrows point to relative maxima for each parameter. 



tion has largely gone untested. Figures 1 and 4 show decreasing dark calcification following 
sunset, reaching zero near midnight followed by an increasing rate of dark calcification 
and an increase in respiration (Figs. 5C and 5D) that rises to a peak at 03:00 well before 
dawn. This pattern has occurred consistently in our mesocosm experiments, with the 
same pattern observed in 30 separate mesocosm runs with different communities under 
various conditions as well as in flume studies at our site (Murillo, Jokiel & Atkinson, 2014). 
Barnes & Crossland (1980) used time-lapse photography to measure diurnal growth in the 
staghorn coral Acropora acuminata and found that night-time extension rate was similar 
to or greater than day-time extension. They suggested that, "symbiotic association permits 
rapid growth because the coral can invest in flimsy scaffolding at night with the certainty 
that bricks and mortar will be available in the morning". Wooldridge (2013) has proposed 
a new model for "dark" coral calcification, whereby O2 -limitation of aerobic respiration 
during the night initiates a homeostatic host response that forms the skeletal organic 
matrix. The matrix formed at night subsequently allows rapid growth of the aragonite 
fibers during the "light-enhanced" period of calcification, when abundant energy derived 
from photosynthesis is available. Perhaps the mid-night calcification minimum observed 
in Figs. 1 and 4 at 00:00 reflects this period of organic matrix formation that precedes the 
03:00 night calcification peak. 
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Figure 5 Diurnal changes in seawater chemistry in the three mesocosms for: (A) total alkalinity Ay, 
(B) aragonite saturation state S2 a rag> (Q pH and (D) dissolved oxygen. 



Diurnal changes in concentration of At, pH, S2 arg and DO 

The variables of At, pH, £2 arg DIC, and DO are concentrations while P ne t and G ne t are flux 
rates. Care must be taken when comparing concentrations to flux rates because flux rate 
can be high when concentration is high or low. Or flux rate can be low when concentration 
is high or low. Figure 4 shows patterns that are difficult to interpret because the figure 
mixes flux rates with concentrations. This issue will be discussed and resolved later in this 
discussion, but first we will compare differences in concentrations of key variables over the 
diurnal cycle (Fig. 5). 

Figure 5 reveals several important patterns: 

1. Alkalinity in the "Algae only" mesocosm remained high during the entire diurnal cycle. 
In contrast, the mesocosms containing corals showed lower Aj (Fig. 5A) caused by 
rapid calcification. Aj reduction by the corals was greatest during the daylight hours 
when G ne t was high (Fig. 1 ) with the difference diminishing during nighttime hours. 
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Figure 6 Hourly DIC versus Aj for the diurnal data measured in each of the three mesocosms. Suc- 
cessive hours are connected by the dotted lines to show the patterns of diurnal hysteresis. 



2. The two mesocosms with algae maintained a higher £2 arag throughout the mid-day 
portion of the diurnal cycle (Fig. 5B) which can be attributed to higher pH resulting 
from rapid rates of algae photosynthesis and coral photosynthesis (Eqs. (6)-(8)), with a 
less pronounced difference during the rest of the cycle. 

3. The extreme difference in £2 arag between the "Corals only" and the "Corals plus Algae" 
mesocosms (Fig. 5B) did not produce a corresponding difference in G ne t between the 
two mesocosms (Fig. 1), which demonstrates that £2 arag is uncoupled from G net and 
explains differences encountered when comparing the £2 arag versus G ne t relationship in 
different systems with different diurnal P net regimes. 

4. Night-time pH and f2 arag values (Figs. 5B and 5C) show less variability than G ne t 
(Fig. 1). 

5. The high biomass in the "Coral plus Algae" mesocosm (Fig. 5D) resulted in the highest 
O2 values during daylight hours (due to photosynthesis) and the lowest O2 during the 
night (due to respiration). The "Algae only" treatment had the second highest daytime 
level O2 due to algal photosynthesis and relatively high levels of O2 at night. 

6. The 03:00 calcification peak observed in Fig. 1 is shown by both a decrease in O2 
concentration and a drop in pH due to accelerated respiration. DO (Fig. 5D) and the 
pH (Fig. 5C) are measured independently and both show this effect to corroborate the 
observation. 

Plotting DIC versus At (Fig. 6) demonstrates the major influence of P net on seawater 
S^arag- Calcification and dissolution shift Aj values horizontally along the abscissa in Fig. 6 
and influence DIC values vertically along the ordinate. However, photosynthesis and 
respiration change DIC along the ordinate without changing Aj. The observed hysteresis 
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pattern results from P net driving G ne t and increasing pH. A linear relationship accounting 
for only half of the variance (R 2 « 0.5) between DIC and At was observed for the two 
rapidly calcifying mesocosms containing corals. This relationship does not hold for the 
low-calcification "Algae only" mesocosm. 

As pointed out by McMahon etal. (2013), connecting the points on a graph of Aj vs. 
DIC reveals a circular hysteresis pattern over the diel cycle as shown for the G ne t versus 
^arag plot (Fig. 6). G ne t can account for changes in both the Aj and DIC concentrations. 
However, P ne t can only account for changes DIC concentration. Therefore, £2 arag is a 
function of the changes in carbonate chemistry due to both P ne t and G ne t> and any 
changes in DIC concentration relative to At will result in different influences on £2 arag . 
For example, in systems with high organic production relative to calcification (Coral plus 
Algae mesocosm), £"2 arag will increase during daylight due to high pH caused by high 
uptake of CO2 used for photosynthesis (Fig. 5D). Conversely, in systems with low organic 
production relative to calcification (Coral only mesocosm), £2 arag will decrease due to the 
uptake of At. Any decrease in Gnet associated with an increase in -Pnet will increase £"2 a raa 
and change the way that G ne t responds to OA. Therefore, any prediction of future global 
changes on coral reef G ne t based on oceanic seawater £2 arag must also take into account 
the influence of future localized P ne t on G ne t and f2 arag as well as changes in carbonate 
dissolution of reef carbonates as described by Murillo, Jokiel & Atkinson (2014) for each reef 
location. 

Comparisons between the hourly and daily G ne t and P ne t values (Table 3) show a 
similar calcification rate for both the "Coral only" mesocosm and the "Coral plus Algae" 
mesocosm in spite of the differences in DO, pH, At (Fig. 5) and P ne t- However, daily P ne t 
for the "Coral plus Algae" mesocosm was only one third of the P ne t of the "Coral only" 
mesocosm. Hourly production was much higher in the "Corals plus Algae" mesocosm 
during the daylight hours, but production was consumed by extremely high respiration 
during nighttime hours. The "Algae only" mesocosm showed very low daily G ne t and 
extremely high daily P ne t- 

Anthony, Kleypas & Gattuso (2011) proposed a model that areas dominated by algal 
beds draw CO2 down and elevate f2 arag , potentially offsetting ocean acidification impacts 
at the local scale. Their model is based on the paradigm that G ne t is controlled by f2 arag . 
Results suggested that a shift from coral to algal abundance under ocean acidification 
can lead to improved conditions for calcification (i.g. increased £2 arag ) in downstream 
habitats and that alga beds can provide a significant mechanism for buffering ocean 
acidification impacts at the scale of habitat to reef. However, this conclusion is at odds with 
the measured values shown in Table 3 and Fig. 5. G ne t in the "Corals plus Algae" treatment 
was the same as the "Coral Only" treatment even though the £2 arag was much higher. In 
addition, the presence of the algae caused a precipitous drop in Pnet* The flaw in their 
model appears to be the assumption that G ne t is controlled by £2 arag . Algal photosynthesis 
increases pH which shifts the equilibrium to higher [CO3 - ] and thus higher £2 arag . There is 
a direct correlation between G ne t and f2 arag for a specific reef community, but not a cause 
and effect relationship. 
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Table 3 Hourly and daily G n et and P ne t values for the three mesocosms. 
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Diurnal changes in material flux (Pnet, Gnet, H + flux and DIC flux) 

G ne t and P net are measures of material flux. DO, pH and DIC are measures of concentra- 
tion. The preceding discussion and numerous publications often compare concentrations 
of one material to flux rate of another material or vice versa. Much more can be learned by 
plotting DIC flux and H + flux rather than [DIC], [H + ] or pH in relation to P ne t and G net . 
DIC flux and H + flux were calculated using the box model and graphed on a 0 to 1 scale 
in the same manner as in Fig. 4 with the result presented as Fig. 7. This figure illustrates 
the dynamic geochemical and physiological relationships involved in coral and coral reef 
metabolism. 

DIC flux (uptake) in the highly calcifying mesocosms containing coral (Fig. 7) increases 
with increasing P net from 06:00 until mid-day peak P net and then decreases rapidly as 
P net decreases with decreasing irradiance. Furla et al. (2000) demonstrated the presence 
of a DIC pool within coral tissues. The size of this pool was dependent on the lighting 
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Figure 7 Plot of normalized data for P ne t, G ne t, inverse DIC flux and H + flux for the experiment. 



conditions, since it increased 39-fold after 3 h of illumination. If we apply this observation 
to the data shown in Fig. 7, it appears that the DIC pool had increased by mid-day, so rate 
of DIC uptake dropped rapidly as irradiance and photosynthesis declined. However, note 
that the high dissipation rates of H + continued for 2-3 h following the peak rates of P ne t 
and G ne t as the corals rid themselves of the backlog of H + generated by rapid calcification. 
Thus the lag of pH behind the peak flux rates of P n et and G net (Figs. 4A and 4B) represents 
a disequilibrium that results from the lag in proton efflux from the corals. The correlation 
between £"2 arag and G ne t is simply the response of the C02-carbonate system to pH as [H + ] 
shifts the equilibria and redistributes the [COj - ] relative to the other DIC components 
of [HCOjT] and [CO2] (Eqs. (3)-(5)) . Therefore £2 arag closely tracks pH whereas G net 
tracks P n et more closely. Changes in f2 arag are a consequence of changes in P net and G net , 
rather than a driver of G ne t- Hence the £2 arag peak and the pH peak lag behind the P net and 
G ne t peaks (Figs. 4A and 4B) due to lag in proton efflux seen in Fig. 7. This observation 
demonstrates the importance of understanding the difference between H + concentration 
andH+ flux. 

During the night the H + flux rate is very responsive to changes in G ne t in the "Algae 
only" and "Coral plus Algae" due to large changes in respiration (Fig. 7). The fluctuations 
of proton flux at night in the "Coral only" mesocosm are dampened considerably 
compared to the "Algae only" treatment. The "Coral plus Algae" mesocosm shows an 
intermediate response. Perhaps the coral skeleton acts as a buffer in a manner similar 
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Figure 8 The flux rate of calcification-dissolution (G ne t) plotted against the concentrations of impor- 
tant variables of the CO2 -carbonate system for the "Corals only" mesocosm. 



to that proposed by Suzuki, Nakamori & Kayanne (1995). The macroalgae lack the large 
skeletal carbonate buffer of reef corals. 

Back to the basics 

The preceding sections have established the importance of using flux rates rather than 
concentrations when we are describing a dynamic metabolic system such as a coral or coral 
reef. Most of the previous research in this area has focused on the relationship between 
G n et> [COj - ] (or its surrogate ^2 arag ), [HCO^~], and [H + ] expressed as pH. Plotting these 
variables in exemplary Fig. 8 is very informative and sheds light on results of previous 
studies. 

A coral must uptake inorganic carbon in order to maintain photosynthesis and 
calcification. As a result [DIC] will decrease no matter which carbonate species (HCOj - , 
CO3 - or CO2) is taken up by the coral (Eqs. (3)— (5)). Thus we see a decline in [DIC] at 
high rates of G net (Fig. 8). [HCO^~], which has been identified as the preferred substrate for 
photosynthesis and calcification (Weis, Smith & Muscatine, 1989; Furla etal, 2000; Roleda, 
Boyd &Hurd, 2012) closely tracks [DIC] during daylight hours. In contrast, [CO3 - ] lags 
behind G ne t and closely tracks pH during the day as shown for £"2 arag in Fig. 4. If [COj - ] 
(or its surrogate ^2 arag ) drives calcification, then how do we explain the lag behind G ne t? 
And if [CO3 - ] is limiting, how do we explain the fact that [CO3 - ] is increasing rather than 
decreasing as the coral calcifies rapidly and takes up inorganic carbon? [CO3 - ] increases 
because of the increase in pH caused by rapid photosynthesis, which shifts the equilibrium 
between [HCO3"] and [CO3"]. Thus, P ne t is the driver of changes in G ne t and [CO3 ] 
(Eqs. (2)-(5)). A basic physiological interpretation of the patterns shown in Fig. 8 is that 
daytime coral metabolism rapidly removes DIC (primarily in the form of HCO^ ) while 
photosynthesis provides the energy that drives G ne t (Fig. 4). Higher pH resulting from 
rapid photosynthesis pushes the equilibria toward higher [CO3 - ]. This scenario results in a 
correlation between G net and £2 arag , with f2 arag as the dependent variable. 
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During the night [HCO3 ], [DIC], [C0 2 3 ~] and pH mirror changes in G ne t- However, 
note that [HCOj - ] diverges from [DIC] and [CO3 - ] diverges from pH in darkness. The 
night divergence can be attributed to respiration causing a decrease in pH. The decreasing 
pH shifts the equilibria so that [CO3 - ] is converted to [HCOj - ], thereby changing the offset 
between the points. This phenomenon is also reflected in the pattern of diurnal hysteresis 
show in Fig. 3. 

CONCLUSIONS AND RECOMMENDATIONS 
Correlations do not establish cause and effect 

Linear regression using £2 arag as the independent variable may be useful as a first 
approximation, but is a poor descriptor of calcification dynamics on coral reefs. Much of 
the existing data on coral calcification was developed in static or low turnover incubation 
experiments under typical laboratory low irradiance artificial light sources on a 12 h light, 
12 h dark cycle (Jokiel, Bohr & Rodgers, in press). This regime results in an unrealistic 
simulation of the actual diurnal cycle that occurs on coral reefs. The standard protocol 
has been to compare linear regressions between or among treatments. Linear regression 
provides a very limited description of the actual relationship between the key factors 
controlling organic and inorganic processes on coral reefs, which are more adequately 
described by data presentations such as that in Table 3 or Fig. 7. The linear regression 
approach does not fully embrace natural diurnal calcification patterns and phase lags 
because these processes are non-linear. The linear regression approach can lead to the 
assumption that £2 arag is the independent variable driving the calcification reaction. 
Nevertheless, correlations at a single location or in a single experiment can result where 
all other factors are held constant because the two quantities are related to some extent. Use 
of £2 arag as an independent variable to compare spatial and temporal variation in G net is 
known to create difficulties (Shamberger etal, 2011; Falter etal, 2012). 

Numerous field and laboratory studies have demonstrated a positive correlation 
between G ne t and f2 arag for corals and coral reefs Erez et al (2011). Well-developed reefs 
occur within a narrow geographic range characterized by open ocean £2 arag > 3.3 (Kleypas 
et al, 1999), which could mean that coral communities have limited capacity to adapt to 
future levels of anthropomorphic ocean acidification (OA) projected for the 21st century. 
Recent reports suggest that healthy coral reefs could cease to exist within this time frame 
as OA continues and oceanic £2 arag decreases (Hoegh-Guldberg et al, 2007; Silverman et 
al, 2009). However, there are inconsistencies in the relationship (slope and x- intercept) 
between G net as a function of £2 arag on various reefs throughout the world (Shamberger 
et al, 2011). For example, Kaneohe Bay, Oahu, Hawaii contains rich coral reefs that 
show extremely high rates of G ne t while living at low £2 arag levels (mean £2 arag = 2.85) 
(Shamberger etal, 2011). Shamberger et al. (2014) report the existence of highly diverse, 
coral-dominated reef communities at the Rock Islands of Palau that are living at low 
saturation states (f2 arag — 1 .9-2.5). These values approach those projected for the tropical 
western Pacific open ocean by 2100 under future OA modeling scenarios. Identification of 
biological and environmental factors that enable these communities to persist at low £2 arag 
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could provide important insights into the future of coral reefs under increasing OA. So how 
do we account for the paradox of rich coral reefs growing at low f2 arag ? Previous work has 
been based on the assumption that G ne t is controlled by or related directly to £2 arag . The 
present investigation indicates that P ne t rather than £2 arag drives G ne t and that calcification 
rate is further limited by proton flux, with pH and Aj playing the major role in controlling 
calcification. f2 arag on coral reefs is simply a dependent variable being controlled largely by 
changes in pH due to photosynthesis. A correlation will exist at a given site, but will not be 
consistent between different coral reef communities. 

The [DIC]:[H + ] ratio correlates with f2 arag in describing G ne t and is useful from a 
physiological point of view because it involves pH and all of the inorganic carbon species 
(Jokiel, 2011a; Jokiel, 2013). However, the [DIC]:[H + ] ratio is simply another variable 
based on concentration (along with pH, £2 arag , CO3 - , etc.) that shows a correlation with 
Gnet- Nevertheless, the [DIC]:[H + ] ratio can be important in describing G ne t in situations 
where £2 arag is decoupled from [H + ] as occurs in the paleo-ocean over time scales greater 
than 10,000 years (Hbnisch et ah, 2012). A calcifying organism must uptake DIC in order 
to continue the calcification reaction and must rid itself of the waste protons. So G ne t 
correlates directly to [DIC] and inversely to [H + ]. 

Results of this investigation (Fig. 7) demonstrate the difficulty that corals encounter in 
shedding waste protons generated during calcification. High rates of H + flux continued 
for several hours following peak G ne t- An important conclusion of this work is that 
measurements of DIC flux and H + flux are far more useful in describing coral metabolism 
dynamics than [DIC] and [H + ] (Fig. 8). Likewise, £2 arag is not a very useful variable in that 
it simply tracks pH (Fig. 4). This pattern becomes clear when one considers that [CO3 - ] 
(and hence f2 arag ) shifts with changing [H + ] as described in Eqs. (2)— (5). DIC flux follows 
Pnet and Gnet and decreases rapidly following peak P ne t and peak G ne t indicating that corals 
can cope more effectively with the problem of DIC supply compared to the problem of 
eliminating H + . 

Future research directions 

The time lag between G ne t and £2 arag reported previously in field studies (Shamberger et 
ah, 2011; Cyronak et ah, 2013; McMahon et ah, 2013) provides evidence that diffusion 
and advection of materials between the coral and the water column involves time delays. 
One reason is that corals convert inorganic carbon to organic carbon, translocate the 
organic carbon to distal calcification sites, store organic carbon as lipid, and can eventually 
convert stored organic carbon back to inorganic carbon (Jokiel, 201 1 b), creating numerous 
possible phase lags for metabolic materials. The second reason for the time lag is that 
rapidly calcifying systems have difficulty dissipating waste protons as shown by continued 
rapid proton efflux for hours after peak calcification (Fig. 7). What other mechanisms 
can account for the phase lag? Boundary layers (BL) can slow the exchange of metabolic 
materials between the coral and the water column. The results of Cyronak et ah (2013) 
revealed that stirring had a net stimulatory effect on A j flux and on the diurnal cycle 
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of hysteresis. Boundary layers slow exchange of metabolic materials, so this is an area of 
investigation that can provide an explanation. 

Three hydrodynamic boundary layers have previously been defined and measured 
(Shashar, Cohen &Loya, 1993; Shashar etal, 1996). The Diffusion Boundary Layer (DBL) 
is only a few mm thick and in contact with the coral epidermis. The Momentum Boundary 
Layer (MBL) controls water movement in the proximity of the sessile organisms and is 
thicker by an order of magnitude than the DBL. The Benthic Boundary Layer (BBL), which 
controls the interactions of the reef with the surrounding sea water, was typically found 
to be more than 1 m thick and characterized by a roughness height of 3 1 cm and a shear 
velocity of 0.42 cm s _1 in the studies. 

The DBL is a thin layer of stagnant seawater adjacent to the coral produced by frictional 
drag. This quiescent layer influences the flux of material between the benthic surface and 
the water column. The transport of Ca 2+ , CO2, CO3 - , HCOj", O2, nutrients and H + 
through the DBL is limited by the physical processes of diffusion and advection (Jokiel, 
1978; Lesser et al, 1994; Kaandorp et al., 2005; Kaandorp, Filatov & Chindapol, 2011). Kuhl 
etal. (1995) found that zooxanthellae photosynthesis resulted in a build-up of O2 in the 
photosynthetic tissue of up to 250% saturation and a tissue pH of up to 8.6 (i.e., 0.7 pH 
units above the pH value of the overlying seawater). In darkness the O2 within the coral 
tissue was depleted by respiration to near anoxic (<2% air saturation) conditions, with 
tissue pH of 7.3-7.4. O2 and pH profiles demonstrated the presence of a 200-300 [xm thick 
DBL that separated the coral tissue from the overlying flowing seawater. Various models 
invoke boundary layer controls on coral metabolism. Kaandorp et al. (2005) and Kaandorp, 
Filatov & Chindapol (2011 ) addressed DBL limitation of DIC influx while Jokiel (2011a), 
Jokiel (2011b) and Jokiel (2013) challenged the paradigm that calcification is limited by 
CO3 - supply on the reactant side of the calcification equation. Rather, he argued that rate 
of dissipation of H + on the product side due to boundary layer conditions can be the actual 
limiting factor. 

Boundary layer limitation of photosynthesis provides an analog to boundary layer 
limitation of calcification. Photosynthetic rate can be limited by rate of waste O2 
dissipation through the boundary layer rather than being limited by supply of reactant 
CO2. By analogy, calcification can be limited by rate of removal of waste protons rather 
than by availability of inorganic carbon. The importance of water motion in reducing 
boundary layer thickness and thereby increasing oxygen flux between the photosynthetic 
organisms and the water column has been demonstrated (Mass et al, 2010). By analogy, 
increased water motion can decrease boundary layer thickness and thereby increase 
removal of protons from the coral. 

Studies of reef metabolism beginning with the classic work of Odum & Odum (1955) 
at Enewetak Reef flat and followed by others (Shamberger et al, 2011; Falter et al, 
2012) were conducted in shallow water reef flats within the BBL in situations where 
unidirectional currents allowed calculation of flux rates. Substantial boundary layers occur 
over all reefs. For example, Price et al. (2012) investigated a range of sites from exposed 
coastal situations to lagoons and found that ambient variability in pH was substantial 
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and oscillated over a diurnal cycle with diel fluctuations in pH exceeding 0.2. Daily pH 
maxima were identified as an important control on calcification. Net accretion among 
sites was positively related to the magnitude and duration of pH above the climatological 
seasonal low, despite myriad other ecological (e.g., local supply, species interactions, etc.) 
and physical oceanographic (e.g., temperature, current magnitude and direction, wave 
strength, latitudinal gradients, etc.) drivers. In general, accretion rates were higher at sites 
that experienced a greater number of hours at high pH values each day. Where daily pH 
within the BBL failed to exceed pelagic climatological seasonal lows, net accretion was 
slower and fleshy, non-calcifying benthic organisms dominated space. Thus, key aspects of 
coral reef ecosystem structure and function are clearly related to natural diurnal variability 
in pH, which is driven primarily by photosynthesis and respiration as P ne t- 

The master variables 

The practice of calculating and comparing linear regressions of G ne t vs. £2 arag to obtain 
a first approximation of calcification rates under different conditions is fraught with 
problems but probably will continue because it is ingrained in science and is convenient 
to use. The correlation of a primary biological response (G ne t) to a primary physical 
chemistry measurement (£2 am g) is attractive, especially in modeling the possible future 
changes on coral reefs. Unfortunately, the physical chemistry concept of f2 arag has 
no basic physiological meaning in describing G ne t other than a correlation with the 
[DIC]: [H + ] ratio (Jokiel, 2013) as well as with other factors such as pH. There is no 
consistent relationship between f2 arag and G net when comparing reefs throughout the 
world (Shamberger et al, 2011). Coral reefs are systems in constant disequilibrium with 
the water column. So we must take care not to be led astray in our thinking about the 
variables that actually drive and control coral and coral reef metabolism and bulk water 
chemistry. The correlation between G ne t and other factors is a result of P ne t driving both 
G net and f^ag (McMahon et al, 2013). The observed phenomenon of diurnal hysteresis 
and diurnal phase lag show the importance of measuring flux rates and emphasizes the 
challenge in predicting the future effects of OA on coral reefs. The method of using linear 
extrapolations of f2 aiag to determine threshold levels that will shift coral reefs from net 
calcifying systems to a net dissolving state has been questioned (McMahon et al, 2013). 
Perhaps predicted changes in ST2 arag in the open ocean can be used to calculate changes on 
reefs if we assume that the baseline on the reefs will change in concert with ocean values 
and that all other processes such as P n et and carbonate dissolution will not be influenced by 
OA. An explanation for the many paradoxes of coral calcification discussed herein has been 
presented as the "Two Compartment Proton Flux Model of Coral Metabolism" (Jokiel, 
201 1 b). This model is focused on localized gradients that influence coral metabolism with 
a focus on proton flux, carbon pools and translocation of fixed carbon. A major feature 
of the model is the presence of boundary layers which control local pH gradients and 
inorganic carbon speciation in addition to proton flux. Results of the present investigation 
support this model. 
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